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Abstract: We present the results of a Bayesian analysis of solar neutrino data in terms 
J> I of z^e — > t'/i.T and Ve Vs oscillations, where Vs is a sterile neutrino. We perform a Rates 



Analysis of the rates of solar neutrino experiments, including the first SNO CC result, and 
spectral data of the CHOOZ experiment, and a Global Analysis that takes into account 
also the Super-Kamiokande day and night electron energy spectra. We show that the 
Bayesian analysis of solar neutrino data does not suffer any problem from the inclusion 
of the numerous bins of the CHOOZ and Super-Kamiokande electron energy spectra and 
allows to reach the same conclusions on the favored type of neutrino transitions and on the 
determination of the most favored values of the oscillation parameters in both the Rates 



and Global Analysis. Our Bayesian analysis shows that Ve I's transitions are strongly 
' disfavored with respect to — > t'^^r transitions. In the case of — > f^^.r oscillations, the 

Large Mixing Angle region is favored by the data (86% probability), the LOW region has 
^ ■ some small chance (13% probability), the Vacuum Oscillation region is almost excluded (1% 

probability) and the Small Mixing Angle region is practically excluded (0.01% probability). 
We calculate also the marginal posterior probability distributions for tan^i? and l^m? in 
the case of v^^r oscillations and we show that the data imply large mixing almost 

with certainty and large values of Am^ are favored (2 x 10~^ eV^ < Am^ < 10~^ eV^ with 
86% probability) . We present also the results of a standard least-squares analysis of solar 
neutrino data and we show that the standard goodness of fit test is not able to reject 
pure Ve — > Vs transitions. The likelihood ratio test, which is insensitive to the number of 
bins of the CHOOZ and Super-Kamiokande energy spectra, allows to reject pure v^, — > Us 
transitions in favor of z/g f^,r transitions only in the Global Analysis. 
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1. Introduction 

The experimental investigation of the solar neutrino problem (see Q) has received recently 
a fundamental contribution from the first data of the SNO solar neutrino experiment Q, 
in which charged-current interactions of solar ®B electron neutrinos with deuterium have 
been observed with a rate about 0.35 times of that predicted by the BP2000 Standard 
Solar Model (SSM) A model independent comparison of SNO and Super-Kamiokande 
Q data show that or Uj- are present in the flux of solar neutrino on Earth at 3.06(7 
^ (see Ref. for a discussion on the probability content of this statement). Since 
the simplest explanation of the presence of or Vj- in the flux of solar neutrino on Earth 
is neutrino oscillations due to neutrino masses and mixing (see [^), this is the second 
model-independent evidence in favor of neutrino mixing, after the up-down asymmetry of 
multi-GeV atmospheric f^-induced events discovered in the Super-Kamiokande experiment 

i- 

In the simplest case of mixing of two neutrinos i^e and i^x that we consider here, we 
have 

f e = cos "dvi + sin 1) V2 , 

Vx = — sm.'d vi + cos-d V2 ■, (1-1) 

where 'd is the mixing angle and z^i, V2 are massive neutrinos with masses mi, m2, respec- 
tively. The flavor neutrino Vx could be either active, x = a with a = ^,t, or sterile, x = s 
(see 0). 

Following the release of the SNO first results, the data of solar neutrino experiments 
(Homestake @, GALLEX SAGE Q, GNO |l|, Super-Kamiokande g, SNO @) 
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have been analyzed in terms of neutrino oscillations in Refs. |13, ^, 14, |l^, |T^, 17, 15, 19 1. 
The above-mentioned 3.06(7 model-independent evidence in favor of the presence of 
or Vr in the flux of solar neutrino on Earth clearly implies that the solution of the solar 
neutrino problem in terms of oscillations of solar i^g's into active u^^ or v-j- is strongly favored 
with respect to transitions into sterile neutrinos Vg- Therefore, assuming the BP2000 SSM 
prediction for the neutrino fluxes produced in the core of the Sun [0] and using a standard 



least-squares analysis, the authors of Refs. ||5|, 14, 15, 1£] calculated the allowed regions in 
the plane of the oscillation parameters tan^t? and Am^ in the case of — ^ i^a transitions 

The standard least-squares method for the analysis of solar neutrino data is an approx- 
imation of the rigorous frequentist method to calculate allowed regions with exact coverage 



|24, 16, 25]. It has been shown in Refs. |24, 25| that an analysis of solar neutrino data 



with a rigorous frequentist method may lead to results that differ significantly from those 
obtained with the standard least-squares method. Unfortunately, the implementation of 
a rigorous frequentist method in the analysis of solar neutrino data is complicated and 
computer-time consuming because it is necessary to calculate the distribution probability 
of the estimator of the parameters for all possible values of the parameters. Moreover, 
several rigorous frequentist methods are available and it is not clear which is the most 



appropriate (if there is one) for the analysis of solar neutrino data (see Refs. [24, 16, 25 1 
and references therein). 

In this paper we present the results of an analysis of solar neutrino data in terms of 
two- neutrino oscillations in the framework of Bayesian Probability Theory (see also [|l^]). 
This analysis has practically the same difficulties of the standard least-squares analysis and 
can therefore be easily performed by the groups specialized in the analysis of solar neutrino 
data. 

The Bayesian analysis of solar neutrino data is not approximate as the least-squares 
analysis in the framework of Frequentist Statistics. Furthermore, Bayesian Probability 
Theory has several advantages over Frequentist Statistics, that have been discussed in 



several books and papers (see |2(;, 27, Eq, 29, pO|). Here we notice only a few facts: 



1. All human statements, including scientific ones, represent knowledge (belief) with 
some degree of uncertainty. Bayesian Probability Theory allows to quantify this 
uncertainty through the natural definition of probability as degree of belief. 

2. Bayesian Probability Theory allows to calculate through Bayes Theorem the improve- 
ment of knowledge as a consequence of experimental measurements. This is how our 
mind works and how science improves. Therefore, Bayesian Probability Theory is 
the natural statistical tool for scientists. 



^The authors of Ref. took also into account the very interesting possibiUty of simultaneous i/e — * Va 
and i^e Vs transitions, that is allowed in four-neutrino mixing schemes (see |2^, ^ ^ and 

references therein). In this case, besides t&v?-d and Am^, there is an additional parameter cos ^ 77 that 
regulates the relative amount of Vc, — > Va and — > Va transitions (0 < 77 < 7r/2). We do not consider here 
this possibility because of its computational difficulties (Bayesian credible regions in the tan^j?-A7n'^ plane 
must be calculated integrating the posterior probability distribution function for the oscillation parameters 
over cos^ rf). 
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3. Probability of an event in Frequentist Statistics is defined as the asymptotic relative 
frequency of occurrence of the event in a large set of experiments. Obviously such set 
is never available in practice. Therefore, Frequentist Statistics is based on imaginary 
data. On the other hand, all inferences in Bayesian Probability Theory are based 
only on the data that actually occurred. 

4. Because of the definition of probability in Frequentist Statistics, the results obtained 
have usually good properties {i.e. they are reliable) only in the case of large data 
sets. In frontier research often a small number of experimental data are available 
(as in the case of solar neutrinos) and the scientist is not interested in long-term 
behavior of inferences, but in getting the best possible inference from the available 
data. Bayesian Probability Theory satisfies this wish. 

5. Results in Frequentist Statistics are based on hypothetical data, i.e. data that could 

have been observed but did not occur. Often different scientists may have different 
ideas on which are the relevant hypothetical data, leading to different conclusions 
(for example, the so-called "optional stopping problem" is well known and discussed 
in the literature). 

6. In Frequentist Statistics there is no way to take into account theoretical and system- 
atic errors, that are not random variables. Nevertheless, since the great majority of 
scientific measurements suffer of systematic errors and scientific inferences depend 
on theoretical errors (they are both crucial in the analysis of solar neutrino data), 
frequentists treat these errors as if they were random variables. Nobody knows the 
meaning of results obtained in this way, with an inconsistent method. Bayesian 
Probability Theory obviously can treat theoretical and systematic errors on the same 
footing as statistical ones, leading to consistent results. 

7. Once a problem is well posed the application of Bayesian Probability Theory is clear 
and straightforward and leads to unique results. On the other hand, in the frame- 
work of Frequentist Statistics several arbitrary and difficult choices that can lead to 
different results are necessary. That is why popular methods (as least-squares) that 
sometimes have poor performances are widely used without a real understanding of 
the motivations. 

8. Since Bayesian Probability Theory allows to calculate the improvement of knowledge 
as a consequence of physical observations, a prior knowledge is necessary. It has been 
argued by advocates of Frequentist Statistics that prior knowledge is subjective, lead- 
ing to undesirable subjectivity in the derivation of scientific results (forgetting the 
often less clear subjective choices of method, estimator, etc. necessary in Frequen- 
tist Statistics). On this problem widely discussed in the literature we want only to 
remark that all human activities, including scientific research, have some degree of 
subjectivity, but communication and collaboration among people working in a field 
allow to reach an agreement on the most reasonable prior knowledge in the field and 
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a way to quantify it. Once the prior knowledge is fixed, Bayesian Probability Theory 
leads to unique conclusions. 



The plan of the paper is as follows. In Section ^ we perform a standard least-squares 
analysis of solar neutrino data. From the calculational point of view this is not an additional 
task, because the calculation of the least-squares function is a necessary step in order to 
derive the Bayesian probability distribution for the parameters. We consider the two cases 
of fe Va transitions (Active)^, and Ve — > transitions (Sterile). In Section ^ we present 
our Bayesian analysis of solar neutrino data. In Section |3.1| we compare the probabilities 
of the models with Va and Ue transitions. In Section |3.2| we present the results 

of our calculation of Bayesian allowed regions (credible regions) in the tan^??-Am^ plane 



for Active transitions. In Section 3.2 we marginalize the posterior Bayesian probability 
distribution for Active transitions in order to calculate the posterior Bayesian probability 
distributions for tan^i? and Am^. Finally, in Section ^ we draw our conclusions. 

In both the standard least-squares and Bayesian analyses we consider first only the 
total rates measured by solar neutrino experiments (Rates Analysis) and then the total 
rates of the Homestake, GALLEX+GNO+SAGE, SNO experiments and the Super-Ka- 
miokande day and night electron energy spectra (Global Analysis). In both the Rates 
Analysis and the Global Analysis we consider also the data of the CHOOZ experiment, 
that are important because they exclude large mixing for Am^ > 10"^ eV^ |3|]. Therefore, 



we consider only values of Am^ smaller than 10 ^ eV^ and we will see that, consistently, 
the allowed regions for the oscillation parameters are limited below this value of Am^. 

In this paper we use the standard names for the regions in ta,v?'&-/S.m?' plane (see ||7|): 
Small Mixing Angle (SMA) for 

10"^ < tan^t? < 10"^ , 10"^ eV^ < Am^ < 10"^ eV^ , (1.2) 

Large Mixing Angle (LMA) for 

10"^ < tan^i? < 10 , 2 X 10"^ eV^ < Am^ < 10"^ eV^ , (1.3) 

LOW for 

10"^ < tan^i? < 10 , 10"^ eV^ < Am^ < 2 x 10"*^ eV^ , (1.4) 
Vacuum Oscillation (VO) for 

10"^ < tan^i? < 10 , 10^" eV^ < Am^ < 10"^ eV^ . (1.5) 



^The treatment of solar Va transitions in the framework of the simple two-neutrino mixing scheme 

in Eq. ( [l.l| ) with Vj: — Va is valid with good approximation in the case of mixing of three neutrinos (vs^ 
i^ii, Vt) with small UeZ as indicated by the results of the CHOOZ long-baseline disappearance 

experiment js^ . In this case Va is a linear combination of and Pt- Furthermore, the treatment is also 
valid in 3-1-1 four-neutrino mixing schemes in which the sterile neutrino is practically decoupled from the 
active ones (1 — |?7s4p <C 1) ^] and U^s is small, as indicated by CHOOZ data p^ . 
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Detection Material and Process 


Data 


37 /"(i. ,, 1 37 ri] . 37 A T. i 
Ul: l^e ~r — ^ Ar + e 

(Homestake g) 


2.56 ± 0.23 SNU 


71Pa- 7/ _i_ 71r"„ , 71p„ _i_ „— 

Lia. z^e + (ja — > KjrG -\- e 
(GALLEX 11 + GNO |12] + SAGE ||]) 


74.7 ± 5.1 SNU 


D2O: Ue + d^p + p + e' 
(SNO i) 


0.347 ± 0.028 


H2O: + e~ — 5- + e~ 
(Super-Kamiokande [^) 


0.459 ±0.017 



Table 1: The rates measured in solar neutrino experiments. The rates of the 

Homestake and GALLEX+SAGE+GNO experiments are expressed in SNU units (1 SNU = 
10"'^^ events atom^^ s^^), whereas the results of the Kamiokande and SNO experiments are ex- 
pressed in terms of the ratio of the experimental rate and the BP2000 Standard Solar Model 
prediction Q. The statistical and systematic uncertainties have been added in quadrature. The 
GALLEX+SAGE+GNO rate is a weighted average of the GALLEX+GNO rate reported in Ref. |l| 
and the SAGE rate reported in Ref. Q . The rate of the SNO experiment is that measured through 
CC weak interactions. 



2. Standard Least-Squares Analysis 

In our standard least-squares analysis (often called "x^ analysis") we consider first in the 



Rates Analysis of Section 2.1 only the total rates measured by solar neutrino experiments 



(see Table ||), and the CHOOZ positron spectra given in Ref. |Q. Then in the Global 



Analysis of Section 2.2 we consider also the Super-Kamiokande data on the day and night 



electron energy spectra. 
2.1 Rates Analysis 

Our least-squares function for the Rates Analysis is written as 



X^ = Xi + X^, (2.1) 



where Xg takes into account the rates measured in solar neutrino experiments and Xq 
takes into account the data of the CHOOZ experiment. 

The X"^ of the solar neutrino rates, Xg, is calculated using the procedure described in 



Refs. [H, H, IJ: 



Ns 

ii j2=i 



fex) 

where Vs is the covariance matrix of experimental and theoretical uncertainties, Rj is 

the event rate measured in the j*^ experiment and R^^^^ is the corresponding theoretical 
event rate, that depends on Am^ and tan^ 6. The indices j,ji,j2 = 1, • • • > Xg with Xg = 4 
indicate the four types of solar neutrino experiments listed in order in Table || together 
with the corresponding event rates and experimental uncertainties. 
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The covariance matrix Vs takes into account the experimental and theoretical uncer- 
tainties added in quadrature^ 32, 36, 24]: 



+ E E (AlnX,)' , (2.3) 

«l,«2=l k=l 

where <t| are the experimental uncertainties given in Table calculated by adding in 
quadrature the statistical and systematic uncertainty for each experiment^. The indices 
hh,i2 = 1, • • • , 8 denote the solar neutrino fluxes produced in the eight solar thermonuclear 
reactions pp, pep, hep, ■^Be, ^B, ^^n, 15q^ i7p^ respectively. The index k = 1, . . . , 12 
labels the twelve input astrophysical parameters in the SSM (5ii, ^33, 5*34, 5'ii4, Sn, 
Luminosity, Z/X, Age, Opacity, Difl'usion, Cr-g^, Sue), whose relative uncertainties AlnX^ 
determine the correlated uncertainties of the neutrino fluxes through the logarithmic 

derivatives 



We adopt the values of an. and AlnXj^ given in Ref. |32|, except for AlnX-j (relative to 
Z/X), whose value has been updated in the BP2000 SSM || from 0.033 to 0.061, and 
Oii2 = Sis, AlnXi2 = 0.181, that have been introduced for the flrst time in the BP2000 
SSM X12 = Sue is the 5-factor for the reaction ^^0(p, 7)^^F that determines the small 
^^F neutrino flux (i = 8). With these values we calculate the fractional uncertainties of 
the SSM neutrino fluxes 



Aln^ssM 



12 



^(a,fcAlnXfe)2 (2.5) 

\ k=l 

listed in Table |2|. One can see that these fractional uncertainties are in good agreement 
with those presented in Table 7 of Ref. |3|. The uncertainty of the hep flux is rather 
small, in view of the fact that no uncertainty is given in Ref. Q because of the difficulty in 
calculating the 5-factor of the hep reaction. However, this is irrelevant for our calculation 
because in any case the contribution of the hep flux to solar neutrino data is negligible (the 
Super-Kamiokande collaboration found that the flux of hep neutrinos on Earth is less 
than 4.3 times the BP2000 SSM prediction at 90% CL, that is one hundred times smaller 
than the *B neutrino flux). 
The quantity 

^(th) ^ ^SSM ^.gh) pi^)(^Am',t^n'e), (2.6) 

such that 

8 

?(th) _ „(th) 



i?f)=5:<\ (2.7) 



^This is the common practice, in spite of the impossibihty to treat theoretical uncertainties in the 
framework of Frequentist Statistics (see item ^in Section ^). 

"^Also this common practice is completely unjustified in the framework of Frequentist Statistics (see 
item 1^ in Section |l]). 
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i 


1 


2 


3 


4 


5 


6 


7 


8 




pp 


pep 


hep 


^Be 




13N 


150 






0.010 


0.016 


0.033 


0.100 


0.176 


0.184 


0.209 


0.239 



Table 2: Relative uncertainties Aln^f^'^ of the neutrino fluxes obtained with Eq. (2^). They 
are in good agreement with those presented in Table 7 of Ref . [|| . 



is the theoretical event rate in the j experiment due to the neutrino flux produced 
in the i*^ thermonuclear reaction in the sun according to the BP2000 Standard Solar Model 
Q. cj^j^^ is the corresponding energy-averaged cross section and p/j^^''(Am^, tan^ 0) is the 
corresponding averaged survival probability of solar z/g's, that depends on Am? and tan^ 6 
(in the case of the Super-Kamiokande experiment, j = 4, also the averaged i^a 
transition probability must be properly taken into account). The quantity AlnCj-**^^ = 

Mjf^'^ /cf^'^ is the relative uncertainty of the energy-averaged cross section cf^^\ The 
values of AlnC^!'*'^^ for j = 1,2,4 (^"^Cl, '^^Ga and H2O experiments) are given in Ref. p2| . 
For the SNO experiment {j = 3) AlnC;^^ ^ = AlnCfg''^ = 3.0 x 10"^ |] and the others are 
zero. 

For the calculation of the probabilities P^^^^Am? ,iaii^ 6) we have used the tables 
of neutrino spectra, solar density and radiochemical detector cross sections available in 
Bahcall's web pages pi, BP2000 Standard Solar Model §. For the calculation of the 



theoretical rate of the SNO experiment we used the charged-current cross section given 
in Refs. |^, The probability of neutrino oscillations is calculated with an unified 

approach that allows to pass continuously from the vacuum oscillation regime to MSW 



transitions |4l|, ^ through the quasi- vacuum regime |^, Q , using the quasi- vacuum 



analytical prescription given in Ref. |4^ (see also Refs. |46, |4^), the usual prescription for 
the MSW survival probability (see [p4| , |4^ ) and the level crossing probability appropriate 
for an exponential density profile [E^, pO] . We calculate the regeneration in the Earth using 



a two-step model of the Earth density profile pl| , 52, 54, 55|, that is known to produce 
results that do not differ appreciably from those obtained with a less approximate density 
profile. 

The of the CHOOZ positron spectra, is calculated as in the analysis A of 



Ref. 1 33 1, with the following approximations. Since we do not know the antineutrino spec- 
trum, the spatial distribution functions of the reactor cores and detector and the detector 
response function linking the real and visible positron energies, for each energy bin we cal- 
culated the oscillation probability at the average energy of the bin and the average distance 



of the detector from each of the two reactors [56]. This approximation is quite good be- 
cause we are interested in values of Am^ below 10~^ eV, for which the energy and distance 
dependence of the survival probability of the i>eS in the CHOOZ experiment is very weak. 
We calculate as in Eq. (13) of Ref. |Q, with the only difference that we neglect the 
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Analysis 


v2 
min 


tan^t? 


Am^ (eV^) 


n.d.f. 


g.o.f. 


Active 
Rates 


8.8 


0.31,3.2 


7.7 X 10"ii 


15 


89% 


Sterile 
Rates 


11.3 


0.36,2.8 


1.1 X 10-1° 


15 


73% 


Active 
Global 


38.4 


0.38 


6.3 X 10"^ 


52 


92% 


Sterile 
Global 


52.8 


0.36,2.8 


4.9 X 10-1° 


52 


44% 



Table 3: Minima X^jj^ of the least-squares function X"^ for Va (Active) or j/g Vs (Sterile) 
in the Rates Analysis of Section ^.1| and in the Global Analysis of Section |2.2| . We give also the 
number of degrees of freedom (n.d.f.), and the corresponding goodness of fit (g.o.f.) assuming for 
^min ^ distribution with n.d.f. degrees of freedom. 



energy-scale calibration factor, whose small uncertainty (1.1%) is practically negligible: 

where Nq = 14 is the number of energy bins and ac is the absolute normalization constant 
with uncertainty a^Q = 2.7 x 10"^ ||3^. We calculate the CHOOZ covariance matrix Vq 
as described in Eq. (12) of Ref. |3^. The only missing information in Ref. [^3| is the value 
of the systematic uncertainties of the positron energy bins, for which only the values for 
positron energy 2 and 6 MeV are given. For the other bins we take systematic uncertainties 
interpolated linearly between these two values. 

The values of the minimum X^j^ of the least-squares function (2.1) and the corre- 
sponding best-fit values of the oscillation parameters for v^, — > Va or fe — > Vs transitions 
are given in the first two rows of Table |. The contribution to Xl^^^ from the CHOOZ 
least-squares function Xq is always 6.0, in reasonable agreement with the minimum value 
(-'^min)cHOOZ = 5.0 and (^no-oscillations) CHOOZ = 5.5 fouud by the CHOOZ collaboration 
in Ref. ||. 

One can see from Table |3| that the minimum X^j^^ lies in the Vacuum Oscillation 
(VO) region for both Active and Sterile transitions. Therefore, there are two minima 
corresponding to the same value of sin^ 2t? (they are connected by log tan^i? — log iaiP'd). 

As written in Table ^, the number of degrees of freedom (n.d.f.) in the Rates Analysis 
is 15, given by 4 solar rates plus 14 CHOOZ bins minus 3 parameters {ac, tan^i?, Am^). 
The goodness of fit obtained assuming for X^j^^ a distribution with n.d.f. degrees of 



freedom (see Refs. [24, ^] for a discussion on the approximate character of this assumption) 
is more than acceptable for both Active and Sterile transitions. 

However, the large value of the goodness of fit in the first two rows in Table ^ looks 
suspicious: if the assumption of a distribution with n.d.f. degrees of freedom for X"^^^ 
is approximately correct, the expected value of X"^^^ is equal to n.d.f. = 15, which would 
give a goodness of fit of about 45%. The small value of X"^ - is mainly due to the inclusion 
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in the fit of the 14 energy bins of the CHOOZ experiment, that contribute with 13 degrees 
of freedom and give a contribution of only 6.0 to -^min- This problem is aheady present in 
the analysis of CHOOZ data performed by the CHOOZ collaboration |33|, where a value 
of (^^in)cHOOZ = 5.0 was found with 11 degrees of freedom. 

From the second row in Table ^ and taking into account that the fit of the four solar 
neutrino rates contributes with 5.3 to the value of ^mirn could conclude that the 
explanation of the solar neutrino problem in terms of fe — > z^s transitions is acceptable. 
This conclusion would be in contradiction with the model-independent evidence in favor of 
the presence of active f^^r in the flux of solar neutrinos on Earth discussed in Refs. ^, |^. 
This apparent contradiction is resolved by noting that the model-independent calculations 
performed in Refs. |2|, ^ ^ are based on the adjustment of the Super-Kamiokande electron 
energy threshold in order to make the response functions of the SNO and Super-Kamio- 
kande experiments to solar neutrinos approximately equal |57, |5^. The authors of Ref. |5| 
found that, given the SNO threshold T^^'^ = 6.75 MeV, the SNO and Super-Kamiokande 
response functions are approximately equal for the Super-Kamiokande threshold = 
8.6 MeV, instead of = 4.5 MeV for which the rate in the fourth row of Table |l] has 
been obtained Q. The low X'^-^^ that we find in the VO region is due to the strong energy 
dependence of the fg Vs transition probability in the case of vacuum oscillations, such 
that the averaged transition probability is quite different in the SNO and Super-Kamio- 
kande experiments if their response function are different. Indeed, for the values of tan^i? 
and Am^ corresponding to X^jj^ we have {Pve^ve)s^o = 0.347 and {Pve^ve)sK = 0.504. 

We have repeated our least-squares analysis using T^^ = 8.6 MeV and the corre- 
sponding rate 0.451 it 0.017 |Q, |5| for the Super-Kamiokande experiment, normalized to 
the BP2000 SSM prediction, and we found Xl^^^ = 14.7 for tan^i? = 1.4 x 10"^ and 
Am^ = 6.0 X 10"^ eV^, i.e. in the SMA region. Now the contribution to ^^jn from the fit 
of the four solar neutrino rates is 8.7, in reasonable agreement with the model- independent 
3.O61T exclusion of pure Vs transitions discussed in Refs. [Q, |5|, ^, which would imply a 
contribution of about 9.4 to ^min- The discrepancy is due to the small remaining difference 
of the SNO and Super-Kamiokande response functions (see Fig. 1 of Ref. [^), such that the 
averaged — > t^s transition probabilities in the SNO and Super-Kamiokande experiments 
are not exactly equal: for the values of tan^t? and Am^ corresponding to ^^j^ we have 
{Pue^ue)'B,^iO = 0.453 and {Pu^^ue)^K = 0.459. Indeed, using "by hand" the same averaged 
probabilities of — > Vs transitions for the SNO and Super-Kamiokande experiments, we 
obtain X^j^^ = 15.5, corresponding to a contribution of 9.5 from the fit of the four solar 
neutrino rates. 

Even using "by hand" the same averaged probabilities of z/g — >■ Vs transitions for the 
SNO and Super-Kamiokande experiments, the goodness of fit is high, 42%. Therefore, 
we see that in the case of Vs transitions the goodness of fit test is misleading, 

because it does not exclude pure Sterile transitions that are excluded at 3.06(T with a 
model- independent comparison of Super-Kamiokande and SNO data @, |5[ 

The goodness of fit test is considered to be appropriate for testing an hypothesis 
against all possible alternatives (see fs^ ). However, in the case of solar neutrino oscillations 
considered here we have two alternative hypotheses: Active Ve — > and Sterile Ve 
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transitions. An appropriate test for alternative hypotheses is the Ukehhood ratio test^, for 
which we find 

Am? ,aQ p(^)(D|tan2??, Am2,ac) ,^ ^, 

Maxta„2^^^^2^„p (D| tan^i?, Am2, ac) 

with the probabihty distribution functions 

p(^)(D!tan2i?, Am2,ac) = ^ \ ^ , , (2.10) 

(27r)^s/2^(2vr)^c/2^' 

where D represents the data, and T = A for v^, — > Va transitions (Active) and T = S for 
fs transitions (Sterile). 
In Eq. ( |2.9| ) we have considered Sterile transitions as the null hypothesis and Active 
transitions as the alternative hypothesis, in order to see if the null hypothesis should be 
rejected or not. This choice is one of the arbitrary choices necessary in the application of 
Frequentist Statistics, whose justification is purely subjective. For example, one may rea- 
son in the framework of 2+2 four-neutrino schemes^ (see [0, 21, 22, ^] and references 



therein) , where transitions of active in sterile neutrinos is predicted in solar or atmospheric 
neutrino experiments. Since no such transitions have been seen in atmospheric neutrino ex- 
periments [30, one could consider reasonable the null hypothesis that the solar neutrino 



problem is due to practically pure i^s transitions''. 

The interpretation of the value of the likelihood ratio for hypothesis testing is another 
subjective choice necessary in the application of Frequentist Statistics, since the likelihood 
does not represent a probability. Usually one chooses an arbitrary value*^, of the order 
of 0.1% or 1%, below which the null hypothesis is rejected in favor of the alternative 
hypothesis. Clearly, the large value of the likelihood ratio (|2.9| ) does not allow to reject the 
null hypothesis of pure Sterile transitions. 

However, we expect better chances to reject the null hypothesis of pure Sterile tran- 
sitions using for the Super-Kamiokande experiment the energy threshold T^^ = 8.6 MeV, 

^ We note, in passing, that in the case of Active transitions the maximum hkelihood is reached in the LMA 
region for tan^i? — 0.30 and Am^ — 1.7x 10~^ eV^. The difference of the location of the maximum likelihood 



and the minimum of (see Table |^ is due to the determinant | Vs | in Eq. ( 2. IC ) , that depends on tan^i9 and 
Am^ . In view of the fact that using the maximum likelihood the best fit of the oscillation parameters in the 



Rates Analysis lies in the same region as the best fit in the Global Analysis (see Section 2.2), a calculation 
of the allowed regions based on the likelihood j2^ may be more robust than the standard calculation based 
on X'^. We do not present it here because we consider the Bayesian analysis discussed in Section ^ definitely 
superior. 

^In order to compare Active and Sterile transitions one must work in the framework of a scheme with 
four or more massive neutrinos. Three neutrinos (ue, Vfj., Vt) are not enough, because the existence of sterile 
neutrinos would be excluded a priori. 

'^A different point of view could be favored by other scientists, as the authors of Ref. [^^ , where a very 
interesting 3+1 four neutrino scheme with practically no transitions of active into sterile neutrinos in both 
solar and atmospheric neutrino experiments has been proposed. In general, in the framework of 3+1 four 
neutrino schemes there may be or not transitions of active into sterile neutrinos in solar and atmospheric 
neutrino experiments. Hence, for solar neutrinos one can consider arbitrarily Active transitions as the null 
hypothesis and Sterile transitions as the alternative hypothesis, or vice versa. 

* Sometimes the value of this small number can be justified in the asymptotic limit of a large number of 
observations (see |59|). Clearly this limit has absolutely no interest in the case of solar neutrino experiments. 
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for which the SNO and Super-Kamiokande response functions are approximately equal Q, 
Indeed, we find 

Maxtan2^,Am^acP^^^(D(^e'' = 8.6 MeV) | tan^^, Am^ ac) _ ^ 
Maxtan2,?,Am2,acP^^KD(rSK = 8.6MeV)|tan2i?,Am2,ac) ' ' 



that is significantly lower than the value in Eq. ( |2.9| ), but still not sufficient to reject the 
null hypothesis of pure Sterile transitions. 

Hence, we see that in the Rates Analysis of solar neutrino data in terms of neutrino 
oscillations Frequentist methods are not able to reject the hypothesis of pure Sterile tran- 
sitions, that are excluded at 3.06a in a model-independent way ||2|, |5|, ^. On the other 
hand, Bayesian Probability Theory allows to assign probabilities to hypotheses and, as we 
will see in Section allows to show clearly that — > Vs transitions are disfavored with 
respect to Ve T^a transitions. 

Let us determine now the allowed values for the oscillation parameters iav?'d and H^m? 
in the case of Ve Va transitions. As done in Refs. js], 15, we do not consider 



T^e — > transitions because they are excluded at S.OGcr in a model-independent way, as 
discussed in Refs. [§, ^, ^. 

In a standard least-squares analysis the 100/3% CL regions in the ta,v?'d-^m? plane 
are given by the condition 

X2<X^j„ + AX2(/3), (2.12) 

where (3 is the Confidence Level (CL) and /S.X'^{(3) is the value of X"^ such that the cumula- 
tive distribution for 2 degrees of freedom is equal to /?: AX2(/?) = 4.61, 5.99, 9.21, 11.83 
for (3 = 0.90(1.64ct), 0.95(1.96o-), 0.99(2.58o-), 0.9973(3o-). 

The allowed regions that we obtain with this standard method in the case of Ve i^a 
transitions are shown in Fig. |l[ One can see that, although the minimum of X"^ lies in 
the VO region, only small areas in the VO region are allowed, because the oscillation 
probability is very sensitive to variations of Am^ and tan^i?. The allowed Large Mixing 
Angle (LMA) region is much larger, because the sensitivity of the oscillation probability 
to variations of Am^ and taiPi} is weak and the value of X^ is low (X^ = 9.1 in the local 
minimum in the LMA region, for tan¥ = 0.34 and Am^ = 2.0 x 10"^ eV^). Also the Small 
Mixing Angle (SMA) and LOW regions are relatively large for the same reason (X"^ = 10.2 
in the local minimum in the SMA region, for tan^t? = 1.4 x 10"'^ and Am? = 8.1 x 10~^ eV'^, 
and X'^ = 12.8 in the local minimum in the LOW region, for tan^i? = 0.60 and Am^ = 
1.4 X 10"^ eV^). The main effect of the inclusion in the analysis of the CHOOZ data is 
to limit the LMA region at large CL below Am^ ~ 10~^ eV^ (in Fig. |^, as well as in the 
following Figs. |, g |). 

2.2 Global Analysis 

In our Global Analysis, instead of the total rate of the Super-Kamiokande experiment 
we consider the data on the Super-Kamiokande day and night electron energy spectra 



presented in Ref. |62|, that contain information on the total rate plus the shape of the 



energy spectrum. The least-squares function is written as in Eq. (|2.lD, with the CHOOZ 
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contribution Xq given in Eq. ( |2.8D . The solar contribution X"^ is written as in Eq. (p.2|), 
but now the indexes ji, j2 run from 1 to A^s = 41. The indexes ji,j2 = 1, 2, 3 refer to the 
first three rates in Table |l|. The indexes ji , j2 = 4, . . . , 22 and the indexes ji , ^2 = 23, . . . , 41 
refer, respectively, to the data on the Super-Kamiokande day and night electron energy 
spectra presented in Table III of Ref . . The covariance matrix Vq, is written as 

11,12 k 
' "j2-3 I l"cor )n V"cor )j2 ' V"cor )n l"cor Jj2 

+ (f^unc^)ii ('^ji(j2-19) +'^ji(j2+19)) > (2-13) 

with 6j = 1 for j > and 9j = otherwise. The quantities for ji, ^2 < 3 are as in Eq. (^). 
For j = 4, ... ,41 we have Uj = {a^^^'^^)j + ((Tunc'')|. The statistical uncertainties {a^^^^^)j 
(for j = 4, . . . , 41) are given in the third and fourth columns of Table III in Ref. l63]. The 
uncorrelated systematic uncertainties (o"uYc'')j (for j = 4, ... ,41) are given by Rj Sj^unc, 
with R^j'^^ listed in the third and fourth columns of Table III in Ref. |62] for j = 4, . . . , 22 
(day spectrum) and j = 23, ...,41 (night spectrum). We assumed that the systematic 
uncertainties of the day and night bins with the same energy are fully correlated. The 
values of (5j,unc are listed in the sixth column of Table III in Ref. ||6^ (we took the bigger 
between the positive and negative values). 

The correlated uncertainties^ (for j = 4, . . . , 41, corresponding to the Super-Kamiokan- 
de energy bins) are divided in theoretical, (iTcor^)ji and experimental systematic, (iTcor 
The theoretical correlated uncertainties (fTcor'')^ = R^^^^ ^fclr come from the calculation of 
the shape of the '^B energy spectrum We list in Table ^ the values of S^^^^^ that 

we have extracted from Fig. 5 of Ref. 1 64 1 . 

The experimental correlated uncertainties are given by {crcm^^)j = R^^^^ ^f'cov Unfortu- 
nately, the relative correlated systematic uncertainties presented by the Super-Kamiokande 
collaboration in the fifth column of Table III in Ref. contains the energy-dependent 
part of the theoretical uncertainty on the neutrino spectrum given in Ref. |64]. In order 
to avoid double-counting, we subtracted quadratically the energy-dependent part of d^^^^j. 

(i.e. — s'^cor^ where j = 4 is the index of the first Super-Kamiokande energy bin) from 
the relative correlated systematic uncertainties in the fifth column of Table III in Ref. p2|] 



(for which we took the bigger between the positive and negative values). The resulting 
values of 6^^^^^. are listed in Table |^. 

This procedure for the treatment of systematic correlated uncertainties is the best 
one that we could find, given the available information. It differs in some way from the 
treatment adopted by other authors |5|, |l4|, 15, 16, 19 1. In any case we would like to remark 
that the effects of these different treatments are very small on the final results, because the 
dominant source of svstematic correl ated uncertainties is due to the theoretical calculation 
^We would like to thank M.C. Gonzalez-Garcia and C. Pena-Garay for explaining to us the treatment 
of correlated uncertainties in Ref. and sharing with us their knowledge on the systematic uncertainties 
published by the Super-Kamiokande collaboration in Ref. that leaded to an improvement of our analysis 
with respect to the first version of the paper sent to the hep-ph electronic archive. 
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of the total neutrino flux, whose uncertainty is about 20% This uncertainty is taken 
into account in the third term on the right-hand side of Eq. ( |2.13D (analogous to the last 
term on the right-hand side of Eq. ( |2.3| )). 

The values of the minimum ^^j^ of the least- 
squares function in the Global Analysis in terms of 
T^e z^a and I'e Vg oscillations are given, respec- 
tively, in the third and fourth row of Table One 
can see that X^jn lies in the LMA region in the case 
of Active transitions, whereas it lies in the VO region 
in the case of Sterile transitions. Again we see that 
the standard goodness of flt test is not able to ex- 
clude pure Sterile transitions. The reason of the low 
values of X^-^^ with respect to the number of degrees 
of freedom is the high number of Super-Kamiokande 
and CHOOZ energy bins whose data do not fluctuate 
enough. This fact could be due to chance or to an 
overestimation of systematic uncertainties. 

However, in the Global Analysis the likelihood 
test allows to reject the null hypothesis of pure Ster- 
ile transitions in favor of the alternative hypothesis of 
pure Active transitions. Indeed, we find 



j 

day 


i 

night 


J, cor 


J, cor 


4 


23 


0.0039 


0.0020 


5 


24 


0.0039 


0.0020 


6 


25 


0.0047 


0.0029 


7 


26 


0.0055 


0.0058 


8 


27 


0.0062 


0.0077 


9 


28 


0.0070 


0.0106 


10 


29 


0.0078 


0.0134 


11 


30 


0.0086 


0.0163 


12 


31 


0.0093 


0.0203 


13 


32 


0.0101 


0.0242 


14 


33 


0.0109 


0.0292 


15 


34 


0.0117 


0.0331 


16 


35 


0.0124 


0.0381 


17 


36 


0.0132 


0.0440 


18 


37 


0.0140 


0.0500 


19 


38 


0.0148 


0.0570 


20 


39 


0.0155 


0.0640 


21 


40 


0.0163 


0.0730 


22 


41 


0.0205 


0.1058 



MaXtan2^,A p^^>{T)\ianH,/\m\ac) 
Maxtan2^Am2,ac p''^K^\ tan^i?, Am2, ac) 



Table 4: 



Values of the relative 



theoretical uncertainties 



^f^L and 
experimental uncertainties S'f^]^^ 



of 



the Super-Kamiokande energy bins 
used in the Global Analysis (see Sec- 
tion O). 



4x10"'', 
(2.14) 

that is small enough. 

Notice that the high number of Super-Kamiokan- 
de and CHOOZ energy bins is irrelevant for the per- 
formance of the likelihood ratio test of alternative hy- 
potheses. In this sense, the likelihood ratio test is 
more robust than the goodness of fit test. 

Assuming that pure Sterile transitions are exclud- 
ed, we present in Fig. ^ the allowed regions in the tan^T?-Am^ plane for Active transitions. 
One can see that there is no allowed SMA region and the allowed VO regions are very 
small. The local minimum in the LOW region has X"^ = 42.1, for tan^i? = 0.64 and 
/S.m? = 1.4 X 10~^eV^, and the local minimum in the VO region has X"^ = 43.9, for 
tan^t? = 0.38,2.6 and Am^ = 4.9 x 10"^° eV^. 

The authors of Refs. [|5, 15, 16| did not find, as us, any allowed SMA region, whereas the 
authors of Refs. |14, 19] found an allowed SMA region at 99.73% CL. These discrepancies 
are due to differences in the treatment of the data and the experimental and theoretical 
uncertainties. Of course, we obtain an allowed SMA region at higher CL, precisely at 
99.98% CL (3.70"). Hence it is clear that the SMA solution of the solar neutrino problem 
is strongly disfavored by present data. We think that further discussions on the precise 
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confidence level at which the SMA region is excluded are not significant in practice and 
their conclusions can change drastically with new data. 



3. Bayesian Analysis 

In this Section we present our results on the fit of solar neutrino data in the framework of 



Bayesian Probability Theory. In Section 3.1 we compare the probabilities of Uf, — > Ua and 



Ve i^s transitions, and in Section 3^ we calculate the Bayesian allowed regions (credible 
regions) for the oscillation parameters Am? and tan^ iD in the case of z/e — Va transitions. 

Bayes Theorem allows to calculate the posterior probability distribution function p(T,- 
tan^t?, Am'^lD, I), where T indicates the type of transitions: T = A for Uf. Ua transitions 
(Active), and T = S for — > transitions (Sterile). Here D represents the data and 
I represent all the background information and assumptions on solar physics, neutrino 
physics, etc. For our analysis Bayes Theorem can be written as 

2n. 2,j.,. p(D|T, tan^^, Am\ I) p{T, tan^^, Am^\l) 
p{T, tan 'd, Am [D, I) = piBlT) ' ^ 

where p(D|T, tan^??, Am^, I) is the likelihood function and p(T, tan^??, Am^ll) is the prior 
probability distribution function (the inclusion of background information and assumptions 
is necessary because in Bayesian Probability Theory, as in real life, all probabilities are con- 
ditional). The probability p(D|I) is known as global likelihood and acts as a normalization 
constant. 

Assuming a normal distribution of statistical and systematic errors, the likelihood 
function is given by 

p(D|T,tan2??,Am^I) = ^ / ^ / dac ^ , , (3.2) 

^' ^ (2vr)^s/2^y '^(27r)^c/2^^' ^ ^ 

where we have marginalized the nuisance parameter aQ. Here A'^g is the number of solar 
data points (A^g = 4 in the Rates Analysis and As = 41 in the Global Analysis) and 
A^c = 14 is the number of CHOOZ data points. Ag is the solar least-squares function and 
1^ is the corresponding covariance matrix, whose calculation in the Rates Analysis and in 



the Global Analysis is explained, respectively, in Sections 2J and 2^. is the CHOOZ 



least-squares function and Vc is the corresponding covariance matrix, whose calculation is 



explained in Section 2.1 



The prior probability distribution function can be written as 

p(T,tan2i?, Am^II) = ^(tan^i?, Am2|T,I)p(T|I) , (3.3) 

where p(T|I) is the prior probability of Vf, — > (T = A) or Vf. Vg transitions (T = S), and 
^(tan^'i?, Am^lT, I) is the prior probability distribution function of the parameters tan^i? 
and Am^ given T and I. 

The prior probability distribution function ^(tan^'!?, Am^jT, I) quantifies the prior 
knowledge on the parameters tan^i? and Am^. Assuming neutrino mixing (included in I), 
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for both T = A and T = S we know that solar neutrino data are sensitive to several different 
order of magnitude of tan^i? and Am^, through vacuum oscillations for Am^ < 10~^ eV^ 
and large mixing angles or resonant MSW transitions for 10~^ eV^ < Am^ < 10^'^ eV^ 
and 10~^ < tan^-!? < 10. Therefore, the most reasonable non-informative prior probabil- 
ity distribution function^'^, that we will use in the following, is a flat distribution in the 
log(tan^'!9)-log(Am^) plane for both T = A and T = S. In this case, using Eq. 
Eq. (|3.1|) becomes 



p(T,tan=*,A™^|D,I)- p(D|T. .an^tf, A»M)p(T|I) 



E 

T=A,S 



dlog(tan2??) dlog(Am2) p(D|T, tan^??, Am^ I) p(T|I) 



(3.4) 

where we have expressed the global likelihood p(D|I) as the appropriate normalization con- 
stant and all probabilities are calculated integrating ^^(T, tan^^?, Am^|D, I) over dlog(tan^t?) 
dlog(Am^). 

3.1 Active Or Sterile? 

In this Section we confront the probabilities of Active Ve — > Va and Sterile Ve — > Vs transi- 
tions using the relation 

p(T|D,I) = y"dlog(tan2i?)dlog(Am2)p(T,tan2t?, Am2|D,I) . (3.5) 

From Eq. ( |3.4D , the ratio of the probabilities of Sterile and Active transitions is given by 
p(S|D,I) _ /dlog(tan2??)dlog(Am2)p(D|S,tan27?,Am2,I) p(S|I) 



p(A|D, I) / dlog(tan2T?) dlog(Am2) p(D|A, tan^i?, Am^, I) p(A|I) 



(3.6) 



Notice that the ratio of the prior probabilities of Sterile and Active transitions factorizes 
out. Since we do not have any prior preference for Sterile or Active transitions, we take 
their prior probabilities as equal, leading to 



j5(S|D,I) _ /dlog(tan2i9)dlog(Am2)p(D|S,tan2T9, Am2,I) 
p(A|D,I) ~ / dlog(tan27?) dlog(Am2) p(D|A, tan^t?, Am2, 1) 

Our result in the Rates Analysis ((D,I) = RatesAnalysis) is 



(3.7) 



p(S I RatesAnalysis) ^ 2 8 x IQ-^ (3 8) 

p( A I RatesAnalysis ) 

It is clear that Sterile transitions are disfavored with respect to Active transitions, in agree- 
ment with the model-independent exclusion at 3.06a" of pure Sterile transitions [Q, 
From the discussion in Section |2.1| , it is clear that the incompatibility of Sterile transitions 
with the data should be enhanced by choosing for the Super-Kamiokande experiment the 



^°We think that this is a prior on which the scientific community can agree, on the basis of the common 
knowledge of the mechanisms of solar neutrino oscillations. Indeed, the authors of Ref. have chosen, 
independently, the same prior. 
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energy threshold Tf^ = 8.6 MeV, for which the SNO and Super-Kamiokande response 
functions are approximately equal ||5|. Indeed, we find 

ff(S|RatesAnalysis,T|K = 8.6MeV) 2 

— — FTfT T = 1.5 X lU . [o.y) 

p(A|RatesAnalysis,r|K = 8.6MeV) ^ ^ 

Notice that, although the number in Eq. ( p.S| ) is comparable with the one in the likelihood 
ratio (2.11), the interpretation is very different. Here we have a ratio of probabilities of 



hypotheses and a very small ratio means that the probability of Sterile transitions is much 
smaller than that of Active transitions. No such interpretation exist for the likelihood 
ratio in the framework of Prequentist Statistics, in which the null hypothesis can be either 
accepted or rejected in favor of an alternative hypothesis on the basis of the comparison 
of the value of the likelihood ratio with a very small fixed number chosen arbitrarily in 
advance (see footnote ^). 

Our result in the Global Analysis ((D,I) = GlobalAnalysis) is even stronger: 

p(S|GlobalAnalysis) ^ ^ ^ ^^.4 
p (A I Global Analysis ) 

that practically excludes Sterile transitions with respect to Active transitions. 

Hence, the Bayesian analysis of solar neutrino data shows that Sterile transitions 
are strongly disfavored with respect to Active transitions, in agreement with the model- 
independent exclusion at 3.06(t of pure Sterile transitions |2|, ^, |6|. We see that Bayesian 
Probability Theory gives an unambiguous and correct result in the comparison of Active 
and Sterile transitions in both the Rates and Global Analysis, contrary to the popular 
goodness of fit test used in traditional analyses of solar neutrino data and the likelihood 
ratio test in the Rates Analysis, as shown in Section ^. 

This better performance of Bayesian Probability Theory in model comparison is first 
of all due to the fact that in Bayesian Probability Theory as in real life one can assign 
probabilities to hypotheses and discuss which hypothesis is more or less favored in compar- 
ison with others, whereas in Prequentist Statistics one is only allowed to either accept or 
reject an hypothesis. Secondly, the judgment of an hypothesis in the framework of Bayesian 
Probability Theory is more robust than in Prequentist Statistics because Bayesian Proba- 
bility Theory allows to estimate the likelihood of a model not only from the best-fit point 
of its parameter space as in Prequentist methods, but from the performance of the model 
averaged over all its parameter space. Thirdly, in the framework of Bayesian Probability 
Theory different hypotheses are fairly compared assigning to them the same prior, whereas 
in the likelihood ratio test it is required to choose a null hypothesis and an alternative 
hypothesis, that are treated quite differently. 

Notice that the Bayesian method does not suffer any problem from the inclusion of the 
numerous bins of the CHOOZ and Super-Kamiokande electron energy spectra, contrary to 
the least-squares method discussed in Section |2[ 

3.2 Credible Regions 

In this Section we present the results of our calculation of Bayesian credible regions (also 
known as highest posterior density regions) in the plane of the oscillation parameters tan^i? 
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and Am^ for Active transitions, that are strongly favored over Sterile transitions, as shown 
in Section 3.1. The credible regions are so-called in order to distinguish them from the 
"confidence regions" calculated in the framework of Frequentist Statistics. Credible regions 
contain a specified fraction of the posterior probability and all values of the parameters 
inside the credible regions have higher probability than those outside. 

The posterior probability distribution function for the oscillation parameters tan^'i? 
and Am^ in the case of Active transitions is given by 

9n . 9,. ^ p(A,tan2??, Am^|D,I) 

tan^^,Am^ A,D,I = ^ n n ' ^.11 

p(A|D,I) 

Assuming Active transitions implies that we take the prior probability of Active transitions 
to be unity: 

p(A|I) = l, (3.12) 
that implies p(S|I) = 0. Therefore, we have 

p(D|I) = p(D, A|I) + p(D, S|I) = p(D| A, I) p{A\l) + p(D|S, I) p{S\l) = p(D|A, I) , (3.13) 

From Eqs. (U), ( |3ll|) , (|]l|) and {^J§ we obtain 

9n . 9,. ^ n(D|A,tan2'i9, Am^,!) 
p{tan^^, Am^lA, D, I) = — ' ' ■ '-^ . (3.15) 

dlog(tan^'(?) dlog(Am^) p(D[A, tan^i?, Am^, I) 

Using the expression (|3.2| ) with T = A for p(D|A, tan^i?, Am^, I), we obtained the 
credible regions with 90%, 95%, 99% and 99.73% posterior probability shown in Fig. ^ for 
the Rates Analysis and in Fig. ^ for the Global Analysis. One can see that these regions 
are similar but larger than the corresponding allowed regions obtained with the standard 
least-squares analysis presented in Figs. |l| and ^. Let us emphasize, however, that the 
meaning of Bayesian and Frequentist regions is quite different and a direct comparison is 
meaningless. 

Bayesian Probability Theory allows to calculate the relative probabilities of the SMA, 
LMA, LOW and VO regions. The posterior probability that the true values of the oscilla- 
tion parameters lie in a region R, with R = SMA, LMA, LOW, VO, is given by 

p(R|A,D,I)= / dlog(tan2i?)dlog(Am2)p(tan2t?, Am2|A,D,I) , (3.16) 
Jr 

where the integration is performed over the appropriate ranges of the parameters given in 

Eqs. (pt)-Cl)- 

In the Rates Analysis we found 

p(LMA|A,D,I) = 0.72, 
p(VO|A,D,I) = 0.12, 
p(SMA|A,D,I) = 0.10, 

p(LOW|A,D,I) = 0.06, (3.17) 
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and in the Global Analysis 



j9(LMA|A,D,I) = 0.86, 
p(LOW|A,D,I) = 0.13, 
p(VO|A,D,I) = 0.01, 

p(SMA|A,D,I) = 10-^. (3.18) 
Hence, the LMA region is favored in the Global Analysis as well as in the Rates Analysis 



(see Ref. |17] for an explanation of the physical reason). This consistency in the Bayesian 



analysis of solar neutrino data is very interesting and promising for future experiments that 



plan to explore the LMA region with terrestrial experiments |65, 66 1. 
3.3 Marginal Bayesian Distributions 



In this Section we marginalize the posterior probability ( 3.15 ) in order to derive the separate 
posterior probability distributions for tan^t? and Am^: 

p{tan^'d\A,'D,I) = J dlog{Am'^) p{tan'^'d, Am'^\A,B,l) , (3.19) 

p{Am'^\A,B,l) = y dlog(tan^i9)p(tan2i?, Am^|A,D,I) . (3.20) 

The posterior probability distributions that we obtained in the Rates Analysis are shown 
in Figs. ^ and ^ and those that we obtained in the Global Analysis are shown in Figs. ^ 
and ^. 

From Fig. |^ one can see that in the Rates Analysis there are two peaks of the posterior 
probability distribution for tan^-!?: one at small mixing angles, tan^t? ~ 10~^, and one for 
large mixing angles, tan^i? ~ 0.3. However, large mixing angles are strongly favored, as 
one can see calculating the corresponding probability: 

p(0.1 < tan^i? < 10|A, D, I) = 0.90 . (3.21) 

In the case of the Global Analysis, Fig. ^ shows that there is only one peak for the posterior 
probability distribution for tan^-i? at large mixing angles, with probability close to unity. 

Figure |6| shows that in the Rates Analysis the posterior probability distribution of 
Am^ has a large and high peak for 2 x 10~^eV^ < Am^ < 10~^eV^, a small and low 
peak for Am? ~ 10"'' eV^, and several sharp peaks for Am^ < 10~^ eV^ (due to the rapid 
oscillations of the transition probability in vacuum as a function of Am^). A calculation 
of the corresponding probabilities shows that large values of Am^ are favored: 

p{2 X 10~^ eV^ < Am^ < 10"^ eV^jA, D, I) = 0.82 , 
^(10^*^ eV^ < Am^ < 2 x 10"^ eV^jA, D, I) = 0.06 , 
p{Am^ < 10'^ eV^|A, D, I) = 0.12 . (3.22) 

From Fig. ^ one can see that in the Global Analysis large values of Am^ are even more 
favored than in the Rates Analysis and very small values are strongly disfavored. Indeed, 
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we obtained the probabilities 

p(2 X 10"^ eV^ < Am^ < 10"^ eV^|A, D, I) = 0.86 , 
p(10"'^ eV^ < Am^ < 2 x 10"*^ eV^|A, D, I) = 0.13 , 
^(Am^ < 10"*^ eV^, D, I) = 0.01 . (3.23) 



4. Conclusions 

In this paper we have presented the results of a Bayesian analysis of solar neutrino data, 
including the recently published rate of CC interactions in the SNO experiment Hj and 



the data of the CHOOZ experiment |33|, that exclude large values of the mixing angle for 
Am^ > 10^^ eV^. We have considered Vf, — > Va oscillations, with a = /i,T, and v^. i/g 



oscillations, where Vg is a sterile neutrino. We have shown in Section 3.1 that the Bayesian 
analysis implies that i^g — ^ i^s transitions are strongly disfavored with respect to fg — ^ J^a 
transitions, in agreement with the model-independent 3.06cj exclusion of pure i^g — ^ J^s 
transitions discussed in Refs. |Q, ^, ^. 

We have also presented, in Section ^, the results of a standard least-squares analysis 
of solar neutrino data and we have shown that the standard goodness of fit test is not able 
to exclude pure Ue — > I's transitions. A problem of the goodness of fit test is the small 
fluctuations of the data relative to the numerous bins of the CHOOZ and Super-Kamio- 
kande electron energy spectra. The likelihood ratio test, that is insensitive to the number 
of bins, allows to reject the null hypothesis of pure Uf, — > transitions in favor of Ug, — > Ua 
transitions only in the Global Analysis. 

In Section we have presented the results of our calculation of Bayesian credible 
regions in the plane of the oscillation parameters tan^-i? and Am? for z^g I'a transitions 
(Figs. ^ and The Bayesian credible regions are significantly larger than the correspond- 
ing least-squares allowed region (Figs. || and ||), presented in Sections 'LI and |2.2| . However, 



it is always necessary to keep in mind that a direct comparison of Bayesian and Frequentist 
regions is meaningless, because they have different properties due to the different definitions 
of probability in Bayesian Probability Theory and Frequentist Statistics. 

Bayesian Probability Theory allows to calculate the probability of separate regions in 
the plane of the oscillation parameters tan^?? and Am^. We found in the Global Analysis 
that the Large Mixing Angle (LMA) region is strongly favored by the data (86% proba- 
bility), the LOW region has some small chance (13% probability), the Vacuum Oscillation 
(VO) region is almost excluded (1% probability) and the Small Mixing Angle (SMA) re- 
gion is practically excluded (0.01% probability). Also the Rates Analysis favors the LMA 
region. 

We have also presented in Section the marginal posterior probability distributions 
for tan^?? and Am?. In the Global Analysis the data imply large mixing almost with 
certainty and large values of Am? are favored (2 x 10~^ eV^ < Am^ < 10~^ eV^ with 86% 
probability). Also in the Rates Analysis large mixing and large values of Am^ are favored. 
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These indications in favor of large mixing and large values for Am^ are very encour- 
aging for future terrestrial experiments as KAMLAND |^5| and BOREXINO [35, 66] that 
have the possibility to explore this region of parameter space. 

In conclusion, we want to emphasize the better performance shown in this paper of 
Bayesian Probability Theory with respect to Frequentist Statistics in the analysis of solar 
neutrino data. In particular, the Bayesian analysis of solar neutrino data is able to disfavor 
clearly Sterile transitions with respect to Active ones, does not suffer any problem from 
the inclusion in the analysis of the numerous bins of the CHOOZ and Super-Kamiokande 
electron energy spectra, allows to reach the same conclusion on the determination of the 
most favored values of the oscillation parameters in both the Rates and Global Analysis. 
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Figure 1: Allowed regions in the standard least-squares Rates Analysis (Section 2.1) of solar 
neutrino rates and CHOOZ data in terms of v^. — > Va transitions. 
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Figure 2: Allowed regions in the standard least-squares Global Analysis (Section 2.2) of the rates 
of the Homestake, GALLEX+GNO+SAGE and SNO experiments, the Super-Kamiokande day and 
night electron energy spectra and CHOOZ data in terms of ly^ — > i^a transitions. 
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Figure 3: Credible regions obtained with the Bayesian Rates Analysis of solar neutrino rates and 

CHOOZ data in terms of Va transitions. The dotted, short-dashed, long-dashed and solid 

curves enclose credible regions with, respectively, 90%, 95%, 99% and 99.73% posterior probability 
to contain the true values of tan^i? and Am^. 
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Figure 4: Credible regions obtained with the Bayesian Global Analysis of the rates of the 
Homestake, GALLEX+GNO+SAGE and SNO experiments, the Super-Kamiokande day and night 
electron energy spectra and CHOOZ data in terms of fa transitions. The dotted, short- 

dashed, long-dashed and solid curves enclose credible regions with, respectively, 90%, 95%, 99% 
and 99.73% posterior probability to contain the true values of tan^^? and Am^. 
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Figure 5: Marginal posterior probability distribution for tan^i? obtained with the Bayesian Rates 
Analysis of solar neutrino rates and CHOOZ data in terms of z/p Va transitions (solid curve). The 
intervals in which the solid curve lies above the long-dashed, short-dashed, dotted, and dash-dotted 
lines have, respectively, 90%, 95%, 99% and 99.73% posterior probability to contain the true value 
of tan^i?. 
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Figure 6: Marginal posterior probability distribution for Am^ obtained with the Bayesian Rates 
Analysis of solar neutrino rates and CHOOZ data in terms of ^ t^a transitions (solid curve). The 
intervals in which the solid curve lies above the long-dashed, short-dashed, dotted, and dash-dotted 
lines have, respectively, 90%, 95%, 99% and 99.73% posterior probability to contain the true value 
of Am^. 
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Figure 7: Marginal posterior probability distribution for tan^i? obtained with the Bayesian 
Global Analysis of the rates of the Homestake, GALLEX+GNO+SAGE and SNO experiments, the 
Super-Kamiokande day and night electron energy spectra and CHOOZ data in terms of — > Va 

transitions (solid curve). The intervals in which the solid curve lies above the long-dashed, short- 
dashed, dotted, and dash-dotted lines have, respectively, 90%, 95%, 99% and 99.73% posterior 
probability to contain the true value of tan^i?. 
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Figure 8: Marginal posterior probability distribution for Am^ obtained with the Bayesian Global 
Analysis of the rates rates of the Homestake, GALLEX+GNO+SAGE and SNO experiments, the 
Super-Kamiokande day and night electron energy spectra and CHOOZ data in terms of — > Va 

transitions (solid curve). The intervals in which the solid curve lies above the long-dashed, short- 
dashed, dotted, and dash-dotted lines have, respectively, 90%, 95%, 99% and 99.73% posterior 
probability to contain the true value of Am^. 
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